Coherent Error Suppression in Multi-Qubit Entangling Gates 
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We demonstrate a simple pulse shaping technique designed to improve the fidelity of spin- 
dependent force operations commonly used to implement entangling gates in trapped-ion systems. 
This extension of the M0lmer-S0rensen gate can theoretically suppress the effects of certain fre- 
quency and timing errors to any desired order and is demonstrated through Walsh modulation of 
a two-qubit entangling gate on trapped atomic ions. The technique is applicable to any system of 
qubits coupled through collective harmonic oscillator modes. 
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The use of spin-dependent forces to create entangled 
quantum systems has become widespread and is 

currently the technique used for the highest fidelity multi- 
qubit operations This powerful technique, first pro- 
posed injd-lsl, has been used to implement quantum algo- 
rithms [9|, create large entangled states 10], test quan- 
tum fundamentals [ill. [l2l. and perform simulations of 
quantum magnetism |l3l. and quantum field theory 
|15| . As these types of experiments are scaled to larger 
numbers of qubits and more complex algorithms, the ac- 
cumulation of gate errors will eventually require quantum 
error correction. Because of the large overhead required 
for quantum error correction, it is important that qubit 
operations be optimized passively in terms of speed and 
robustness to non-ideal control environments. In this pa- 
per, we show how ideas similar to the spin-echo pulse 
sequence [3] and those developed in the context of dy- 
namical decoupling [l7l - [l9j can be used to optimize the 
M0lmer-S0rensen (MS) gate that is based on the spin- 
dependent force. 

Spin-dependent force gates operate by coupling the 
qubit states to a collective external degree of freedom re- 
ferred to as a quantum bus. The coupling is switched on 
for an amount of time that introduces a particular phase 
between the spin states and at the same time leaves them 
disentangled from the external degree of freedom at the 
end of the gate. While the relative spin phase is fairly 
robust due to geometric features [20|, the disentangle- 
ment of the qubit space and the quantum bus at the 
end of the operation may be more susceptible to exper- 
imental errors and is equally crucial to achieving a high 
fidelity gate. Imperfect timing caused by noise on the 
energy splitting of the qubit can be suppressed by the 
insertion of an additional swapping pulse on the qubit 
states in the middle of a two qubit gate [21|] or, as pro- 
posed in j22], by a TT phase shift in the drive field. In 
this Letter, we present new analytic results which gen- 
eralize these ideas and show how frequency and timing 
errors can theoretically be suppressed to any desired or- 
der with an optimized gate sequence that does not rely 
on the insertion of additional tt pulses within the gate. 



Furthermore, the technique is demonstrated using atomic 
hyperfine qubits driven by a stimulated Raman process 
and shown to perform much better than the standard op- 
eration described in the original proposal [ij. Similar to 
the single-qubit composite pulses |,23j originally designed 
for error suppression in NMR experiments now being in 
widespread use in other quantum information systems, 
this composite pulse should be applicable to any system 
of qubits coupled to a driven harmonic oscillator such as 
superconducting flux qubits 24] or cavity QED 22 [. 

In trapped-ion systems, the spin-dependent force cou- 
ples internal atomic states of neighboring ions through 
the collective modes of motion generated by the Coulomb 
interaction. In the MS scheme, a spin-dependent force 
is created by off-resonantly driving the flrst-order red 
and blue sideband transitions simultaneously. The in- 
teraction Hamiltonian takes the form H — il/2{a+e^^' + 
(T_e~''^0(ae"'**e*"*'" + ate'^'e"''^") where ft is the side- 
band transition frequency, a± are the raising and low- 
ering operators for the qubit, |a^,a} are the creation 
and annihilation operators for the collective harmonic 
oscillator mode, and S/2Tr is the symmetric detuning of 
the drive field from the sidebands [6]. The sum phase 
fps — (0b + 4'r)/'2 of the red and blue sideband drive 
fields determines the eigenstates of the spin operator in 
H. The difference phase 0™ = (c/ib — 4>r)/'^ determines 
the phase of the time-dependent displacement of the mo- 
tional state. For the general case of N ions, the time- 
evolution operator is given by, 

U{t) = f,-^ lo dt'H{t')-^ dt' S;^' dt"[Hit'),H(t")] f^^^ 
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spm operator is given hy = 
'''+'* 6**^° + (J^l^ e~'''^% the time-dependent displace- 
ment coefficient is a{t) ^ n/2 di'e"*('')*'e*'^" and $(t) 
is a time-dependent phase that depends only on and S. 
When a collection of trapped ions that are each identi- 
cally prepared in an eigenstate of az evolves according to 
([2|), the spin-dependent displacement operator splits the 
motional wavepacket into TV -1-1 pieces that execute circu- 
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lar trajectories in phase space according to the definition 
of a{t). The term in ([2]) proportional to 5*^ imprints 
a relative phase on the various spin states, allowing the 
operation to be used as an entangling operation. 

In order to prepare a pure spin state with this type of 
operation, the entanglement between the spin and motion 
must disappear at the end of the gate. When the gate 
time tg is not equal to 27rj'/i5, where j is any non-zero inte- 
ger, the motional wavepackets do not trace out closed tra- 
jectories in phase space and therefore will not be fully dis- 
entangled from the spin state. The required level of preci- 
sion grows with higher temperatures since the overlap be- 
tween two states separated in phase space decreases expo- 
nentially with temperature. To see this, consider a qubit 
under the influence of the time evolution operator in ([2|) . 
If the initial motional state is assumed to be a Gaus- 
sian state ipi^) with an uncertainty in position Ax and 
we describe a small timing or detuning error in the gate 
operation as an unintentional momentum displacement 
hq, then the overlap between the two motional states is 
given by J^^dxip*{x)e~'"''^ij}{x) = exp[-^{qAx)'^]. For 
a harmonic oscillator in a thermal state, Ax increases 
approximately as Vt for ksT > hw meaning that the 
overlap between the two states decreases exponentially. 
As shown in Fig. [ija), a small detuning error can largely 
be corrected with a second pulse whose phase has been 
shifted by tt. We now discuss how to generalize this sim- 
ple pulse sequence in order to suppress larger errors of 
this type. 

Suppose there is a symmetric error A in the detuning 
such that 5 = 2'K/tg + A that could be the result of a 
change in the trapping frequency. The error in the oper- 
ation results in some residual entanglement between the 
spin and motion that can be quantified by the magnitude 
of ao{tg) = n/2 /q" dte-'^^+^^\ (which goes to zero for 
A = at tg = 2-Kj /S). We will show that by switching 
either 0s or <prn between and tt at times prescribed by 
certain Walsh functions, the effect of A on the magni- 
tude of a{tg) can be suppressed to any order. A Walsh 
function, denoted here as W(fc, a:), is a piecewise con- 
stant function that alternates between the values ±1 at 
certain values of x depending on the dyadic-ordered in- 
dex k [i^, (see Fig. mb)). If 4>r and ^f, shift together 
between and tt, then 0^ shifts between and tt, but 0^ 
remains constant and can be assumed to be without 
the loss of generality. Note the effect of the phase shift 
= 0s = TT is equivalent to shifting the motional 
phase 4>m = => 0„i = while keeping 0s constant. 
Both of these phase shifts are equivalent to the mapping 
H ^ —H, which can also be achieved with tt pulses on 
the qubit states as done in [21|. Although the tt phase 
shifts and tt rotations are ideally equivalent, the phase 
shift switching time and precision is limited by electron- 
ics whereas the microwave rotations depend on qubit con- 
trol that might be subject to the same noise source that 
generates A. If the times at which phase shifts occur is 
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FIG. 1. (a)The phase space trajectory of the motional 
wavepackets of a single ion during a W(l, t/tg) spin-dependent 
force operation with a small detuning error, A/(5 <gi 1. The 
solid and dashed curves show the two different trajectories 
taken by the two different wavepackets associated with spin 
up and spin down in the Si basis with only the spin up tra- 
jectory being labeled for clarity. After being initialized to 
a state centered at the origin, the spin-up wave-packet be- 
gins its clockwise motion near the point labeled 1. Halfway 
through the operation, near the point 2, the phase of the drive 
field is advanced by tt, changing the direction of the applied 
force. At the end of the gate, near point 3, the wave-packet 
ends up much closer to the origin than the turning point near 
point 2. Therefore, the two wave-packets have more over- 
lap at the end of the gate which is the key to achieving a 
higher fidelity operation. (b)The Walsh functions W(l,a;), 
W (3, x) and W (7, x) are shown. Notice that W(3, x) can be 
constructed as two sequential W(l, x) functions with a phase 
fiip on the second pulse. Likewise, W(7, x) can be constructed 
as sequential W(3, x) pulses with a phase fiip on the second 
pulse. 



determined by the zero crossing times of W (fc, t/tg), then 

Sn = '^{k,t/tg)J2f <J+ + <J- = W{k,t/tg)XN- When 
modulating 0s in this manner, the new displacement op- 
erator in ([2]) is given by, 

(3) 

By choosing a Walsh function with index k — 2" — 1 where 
n is an integer and a detuning S — 2"+^7r/ig, phase flips 
only occur at integer multiples of 2tt/6 and the effect of 
A can be suppressed to any order. This statement rests 
on the following equality. 



.1 n 

/ da;W(2" - 1, 2;)e±*2"+V:r ^ ^^^i 

•'0 1=0 



(4) 



where a/ is a constant, (proof in the supplemental mate- 
rial). If the function e^*'^' is expanded in a Taylor series, 
the identity in ([4]) ensures that the displacement opera- 
tor wiU be given by D2"-iitg) = i + 0(A"+i), where I 
is the identity operator. 
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FIG. 2. A single ion prepared in |^) is subjected to the stan- 
dard and composite spin-dependent force operations and then 
measured in the o"z basis. The data shown are plotted to- 
gether with theoretical curves assuming an initial thermal 
state of motion with the average excitation number n = 7. 
(a) The data show the probability of finding the ion in |t) as 
a function of the symmetric detuning 5 for tg — 100 usee. On 
resonance, 5 = 0, the motional wave-packets quickly become 
entangled with the spin state, resulting in a maximally mixed 
spin state. For finite 5, the wavepackets trace out circles in 
phase space resulting in partial revivals of the initial spin state 
when 5tg/2n is a non-zero integer, (b) The spin-dependent 
force operation is implemented using W(l,t/V2tg) for the 
phase (l>s- (c) W(3, t/2tg) is used for the phase (ps- Note the 
narrow resonance at Stg/2Ti — 2 corresponds to a trajectory 
where the phase flips occur when the motional wavepackets 
are not at the origin. 
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FIG. 3. The state fidelity of a two ion MS gate as a function of 
the detuning 6 is compared for the first three Walsh functions 
being used for (j>s- In all three data sets, the ions were side- 
band cooled to the motional ground state before implementing 
the gate and tg = 0.3ms. In figures (a), (b) and (c) the mea- 
sured state fidelity as a function of the detuning 5 is compared 
with the theoretical curves. To account for other imperfec- 
tions in the experiment, the theoretical curves have an overall 
scale factor. The data shows a maximum fidelity of ~ 0.88 for 
the standard pulse, ~ 0.90 for the W{l,t/\/2tg) sequence and 
~ 0.93 for 'W{3,t/2tg). The characteristic width of the high 
fidelity region is also clearly larger for the higher order pulse 
sequences as predicted by the theory, (d) shows the fidelity 
curves for W(0)(blue), W(l)(red), W(3)(green), W(7)(purple) 
and W(15)(cyan). The curves are shifted in frequency in order 
to facilitate comparison and highlight the increasingly large 
regions of high fidelity for the higher order sequences. 



To demonstrate the power of the composite pulse se- 
quence, we use a qubit defined as the clock states in 
the S'i/2 hyperfine manifold of a Yb"*" ion in an RF 
Paul trap which can be initialized and read out us- 
ing the techniques described in [2^. These states, 
{|F = 0, mF = 0) = li), |F = l,mF = 0) = It)}, have a 
splitting of 12.6428 GHz and can be coupled to each other 
directly or through the harmonic oscillator using stimu- 
lated Raman transitions. As described in 27[, the Raman 
transition induced spin-dependent forces are created by 
the beat notes between two optical frequency combs that 
are generated by a 355nm mode-locked pulsed laser. The 
UV pulses have a duration of ~ 10 psec at a repetition 
rate of 80.57 MHz and have a center wavelength that is 
optimal for minimizing off-resonant scattering from the 
excited P states in Yb"^ [sO] ■ At the position of the ion, 
the two beams are cross-polarized and mutually orthogo- 
nal to a magnetic field of 5 G with a geometry such that 
the momentum kicks associated with the stimulated Ra- 
man process only excite the transverse modes of motion 
which have a resonance frequency of 1.5 MHz. The two 
Raman beams are frequency shifted with AOMs to set 
up the appropriate beat notes in the interference field 
at the location of the ions. As described in 27[, driv- 
ing one AOM with a single frequency and the other with 



two frequencies generates the bichromatic beat note that 
gives rise to the MS interaction. The red and blue phases 
(j)r/b are therefore defined by the phases of these two RF 
drive frequencies. The composite pulse is implemented 
by splitting the operation into segments, between which 
the phases (j)r and (j)b are shifted. In this setup, symmet- 
ric detuning errors might be the result of fluctuating RF 
trap voltages which manifests itself as noise on the oscil- 
lation frequency. Asymmetric detuning errors arise from 
a change in the qubit splitting which might stem from a 
noisy local magnetic field. 

The effect of Walsh modulation on the spin-dependent 
force can be plainly seen with a single ion. In the case 
of a single ion, the phase is global and the only 
operation that results in a pure spin state is one that 
simply restores the initial spin state. The disentangle- 
ment of the spin and motion and consequential revival 
of the initial spin state should occur when 5tg /2tt = 2"j 
where 2" — 1 = /c is the Walsh function index. If the 
spin is initialized to \\.) in the Gz basis and the joint 
spin-motion state after the operation is p, then the state- 
fidelity is Fi = Tr [14,) (tip]. Ignoring heating effects and 
assuming an initial thermal state of motion, the fidelity 



IS 



1 + exp 



-(n + 1/2) \2au{tg)Y 



where n is 
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the average excitation number of the harmonic oscihator 
and ak{tg) = § /g" dm {k,t/tg) e''^^+^'>K Because the 

lowest order term for the infidelity is Odafej^), the infi- 
delity of the Walsh modulated operation at 5tg/2-K = 2"j 
is 0(A^"+^). This effect is clearly seen in Fig. [2] where 
the higher order Walsh sequences exhibit spin revivals of 
high purity over a much larger range of detunings. The 
increase in the width of the spin revival regions for the 
higher order Walsh sequences is evidence of the higher 
tolerance of small detuning errors. Note the same argu- 
ments apply to the case of small timing errors and errors 
resulting from an asymmetric detuning from the red and 
blue sidebands such as a change in the qubit splitting. 

The effect of the Walsh modulation on a two-qubit gate 
is more complicated than that of a single qubit operation 
since the term in ^ proportional to S*^ must be taken 
into account. The Walsh modulation of 0^ changes the 
evolution of in general, but not in the case where 
5 = 2"+^7r/tg since the evolution is a series of closed cir- 
cles in phase space. In this case, = and a 
fully entangling operation is achieved when $(tg) — tt/2. 
This implies that in order to use W {2^ — l,t/tg), the 
gate time must be at least tg = 2"/^7r/f2. While the ex- 
ponential nature of this composite gate becomes daunt- 
ing for large n, small errors can easily be corrected 
with a modest increase in the gate time. In the case 
of two ions, the maximally entangling operation ideally 
implements the transformation |||) |||) + e'^|tt) 
where the phase 9 is determined by the phase of the 
drive field. With this target state, the fidelity is F2 — 



1/4 



3-(fi+l/2)|2Qfc(t,)|= _|_ ^g-in=*fc(tg) 



and is measured 
in the same manner as described in The phase 

Mtg) = E->,=o Im y*^,] -^{tg~l Eto MSU)) is 
written here in terms of sums over the different parts of a 
pulse sequence. The parameters tt refer to the duration 
of the {i + segment of a sequence and the parame- 
ters if, = (-)* J^'_^ dte-'^* with t_i = and tu = tg. The 
data shown in Fig. [3] compare the state fidelity of two 
ions using the composite pulse sequences for the maxi- 
mally entangling gate and shows that the standard pulse 
sequence is outperformed by the Walsh modulated pulse 
sequences in terms of both the maximum fidelity and the 
characteristic width of the high fidelity region. The value 
of r2 is used as a fit parameter for the theoretical curves 
with the same value, 27r x 1.47 kHz, used in all three 
plots. 

Walsh functions have long been known by the elec- 
trical engineering, astronomy and radio communications 
communities to have useful error correcting properties 
(25} . While the Walsh functions are not the only op- 
tion for choosing how to modulate the drive field of the 
spin-dependent force gate, we hope their introduction in 
the context of quantum control provides a useful tool for 
the further development of dynamical decoupling and re- 



lated areas. In the formalism of dynamical decoupling, 
the function ak(tg) can be viewed as an optimized fil- 
ter function designed to suppress the effects of a noise 
source centered at 6/2^ [28|. Walsh modulation is opti- 
mal in a different sense than that of the Uhrig sequence 
which can, with an increasing number of pulses, sup- 
presses higher moments of the noise spectrum in the zero 
frequency limit 18|. While the Uhrig filter is optimal 



in the number of pulses used for a given order of noise 
suppression, the Walsh filters are optimal in the number 
of elementary sequences, (see supplemental material for 
the definition of an elementary sequence). This minimal 
number of elementary sequences not only allows for a 
simple mathematical construction, but also means that 
these functions are easy to synthesize using simple inte- 
grated circuits. 

By introducing the idea of a Walsh modulated spin- 
dependent force, we have shown theoretically and ex- 
perimentally that it is possible to suppress errors that 
are linked to the residual entanglement between the spin 
and motion, thereby alleviating the required precision of 
the control fields. Because the detrimental effect of any 
error increases exponentially with the initial tempera- 
ture of the harmonic oscillator, the technique may also 
decrease the amount of resource intensive cooling that 
must be done in order to achieve a high fidelity opera- 
tion. As quantum information experiments progress, this 
technique of coherent error suppression in quantum bus 
operations might prove to be an important ingredient in 
scaling toward larger systems and more complex algo- 
rithms. 
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Supplemental Material 



In this section, we aim to prove the following identity 
by induction, 



„1 n 

•^0 1=0 



x' = 0. 



(5) 



The construction of the Walsh functions is simple in 
terms of the elementary sequences known as Rademacher 
functions R(n,i) = sign [sin (2"7ri)]. The dyadic order- 
ing of the Walsh functions allow them to be defined 
in terms of the Rademacher functions as W (n, t) = 
nS R(»:i) ^hen n is expressed as a binary num- 
ber n = &™2'" + .... + 6o2° and 6^ = or 1. With this 
definition, it is easy to see that choosing the index 2" — 1 
for the Walsh function means that all the hi coefficients 
are 1. We now prove the base case, n = 1. 



nl 1 

l-i 1 

/ dxR(l,a;)e*'*'^^^a,x' 

1=0 

Jo 

.1/2 

ai / dxe'^'''-' {x^ {x + l/2))^{). 
Jo 



For the inductive step, assume that ^ is true and look 
at the 2"+i - 1 case: 



r 1— 

/ da;W(2"+i - 1, x)e'^"*''"' V 
"'0 , n 



n-l-l 

+ 2^^ \ ^ I 

aix 



n+l 



= / dxR{l,x)....R{n + l,x)e'^" ^ a,x' 
■^0 1=0 

, „2 n+l n+l J 

= 1 dxY[R{^,x/2)e^'"''-^J2-^{?) 

•^0 i=l 1=0 

„i n+l n+l J 

= 1 dxl[R{^,x/2)e^'"''-^Y.-'{^) 

■^0 i=l 1=0 

+ 1 dxY[R{^,x/2)e^'"''^^J2'''{^) ■ 

i=l 1=0 



For the next step, note that for an integer ri > 0, 
R{n + l,x/2) = R(n,x), which allows the expression to 
be written as. 



2 / "^^j J|R(i,a;)e 

2 / '^2: J|R(i,x)( 

o / '^a:]TR(i,a;)« 
"^•^1 »=i 



n+l 



i2" + Va; 



n+l , 



i2" + Va; 



,i2"+Vx 



i=0 

n+l 

E 

n+l 

E 

1=0 



ai 



In the next step, the substitution x' = x — 1 is made and 
we take advantage of the fact that R(n, x + 1) = R(n, x) 
for n > 1. 



1 



n+l 



da;W(2"-l,a;)e^2"+Vx^^^i 



rfa;W(2" - l,x)e 



dxW{T - l,x)e' 



2' 



(=0 

n+l ; 



E^Eu'^'^ 

1=0 k=0 



m 



/=0 



which is zero by assumption since 6; is a constant, thus 
concluding the proof. If the Radamacher functions are 
considered to be elementary sequences, then the Walsh 
filter is optimized in this resource for the task of sup- 
pressing the errors in the spin-dependent force operation 
that are discussed in this paper. 



